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In this paper, shunt active power filter (SAPF) is designed to address the 
problem of current harmonics in the source current arising from a nonlinear 
load and improve the quality of electric power. That will be by compensating 
reactive power and harmonic currents. The PI controller responsible for DC- 
link energy storage tuning was developed using the Lévy flight distribution 
algorithm (LFA). It is a novel, previously unused optimization approach to 
suggest relative gain Kp and integration gain Ki gain values for the PI 
controller. This approach aims to get the best dynamic performance of SAPF, 
speed up the convergence rate, get the fastest best stability and bypass constant 
voltage advancement for DC-link. The model was tested and implemented in 
MATLAB simulation software. The result of total harmonic distortion THD 


p-q theory showed the efficiency of this method compared with the results of traditional 

PWM PI. The design of the model in the current reference frame was based on 

SAPF instantaneous reactive power (pq) theory. 
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1. INTRODUCTION 

The increasing use of non-linear loads represented by electronic transformers and various electronic 
devices led to major problems in the power system, corresponding to the increase in harmonics. It negatively 
affects the quality of electrical energy, a decrease in the power factor, and poor efficiency of electric power 
[1]. Harmonics of the higher order of the power lead to additional losses such as heat problems, turbulence, 
and equipment failure [2]. A reason for the formation of harmonics in non-linear loads is represented in 
electronic power transformers. As a result of the properties of these electronics in crystal vibration and the 
effect of the fast-switching process which arises from this fast-switching angle (a). As a result of the dynamic 
and electromagnetic fasting response in each oscillation through the wave during the operation process, due to 
the greatest and increasing need for power electronics devices [3], it is flexible in operation and less 
maintenance. Therefore, it is necessary to take procedures, methods, and a mechanism for dealing with the 
harmonics resulting from these devices. The passive power filters can be used to remove harmonics of the 
electric current wave of the source resulting from these devices and solve their problems; this method has many 
negatives and problems such as resonance, electromagnetic echo, increased reactive power, and others, to avoid 
the problems associated with the use of passive filters, the use of active power filters (APF) were used, which 
represents the best solution for treating harmonics in the electric current, and it is better than traditional passive 
filters in this performance [4]. The shunt active power filter (SAPF) helps in correcting the current waveform, 
reducing harmonics, obtaining a wave as close as possible to the sinusoidal shape, compensating the reactive 
power, and eliminating unwanted frequencies, and therefore; its use has increased widely [5]. The SAPF device 
is connected in parallel with the network between the source and the load at the common coupling point (PCC), 
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as shown in Figure | and works on the principle of harmonic sensing in the source current wave and generates 
a compensated current injected into the network to get the best results we use different control strategies as 
indicated by the literary references and their development Including methods for calculating power, current, 
and pulse gate generation techniques for pulse width modulation (PWM) of the DC-AC inverter, to achieve 
the proper compensated harmonic current (SAPF) injection into the network, it is necessary to track the source 
current and filter current and maintain the voltage via DC-link taking into account the use of time-domain as 
well as frequency domain. The compensated current is generated by analyzing the distorted voltage and the 
source current signal and separating the harmonics using real-time analytics [6]. A time-domain approach is 
used for ease of use in the design of the controller compared to the frequency domain. The current reference 
frame is based on instantaneous reactive power (pq) theory, which is one of the control methods in the time 
domain approach, where the voltage difference and harmonic current are calculated to estimate the 
compensated reference currently. Used for power and current calculations. It is one of the commonly used 
methods. The strategy of controlling source voltage (DC) of the inverter and comparison with constant voltage 
is using PI, this model was developed by optimizing the performance of the PI controller, which is responsible 
for tuning the energy storage in DC-link, using the Lévy flight distribution algorithm (LFA), a new optimization 
method for suggesting Kp and Ki gain values is suitable for the controller to improve the dynamic performance 
of the SAPF and maintain the required voltage level via DC-link and to obtain the best possible result for total 
harmonic distortion THD compared to the traditional PI control system and various improvement methods for 
other researchers. The model was tested and the results verified in MATLAB simulation software. The system 
under study is a three-phase three-wire system consisting of a three-phase voltage source for the alternating 
current and non-linear load that represented by electronic devices or electronic power inverters, shunt active 
power filter (SAPF), and network impedance as shown in Figure 1 [7], [8]. Active power filter consists of a 
PWM inverter type VSI voltage source inverter fed from a DC source. The APF is connected in parallel 
between an alternating voltage source and a load at the common coupling point (PCC) in which the equation 
of the following currents is fulfilled: 


=l -liş (1) 


where I; represent filter's current, I, is load's current and Is is source's current. 
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Figure 1. Block diagram of the SAPF 


2. CURRENT REFERENCE GENERATION TECHNIQUES 

In this study, the active power filter (SAPF) was designed by adopting the pq reactive power theory 
(Akagi, Kanazawa, and Nabae in 1983-1984) in power calculations and generating reference currents [9], [10]. 
This method is valid for the operation of the SAPF device in the transitional or steady-state and general voltage 
state because it depends on the time-domain approach and simultaneous detection method as shown in Figure 2. 
Allow the SAPF device to control in real-time because of the simplicity of the calculations [11]. This method 
includes only algebraic equations. The compensated harmonic current is generated by analyzing the buffer voltage 
and current distorting harmonics and separating the harmonics using instantaneous transforms [12]. 

The calculations of the pq method include only algebraic equations for the counting of the 
compensated reference current arising in a SAPF device built in a three-phase system. Below, shows the 
equations of Clark transformations for voltage source and load current in coordinates a b c to binary coordinates 
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in the following relations [13]. Clark's calculations method uses the two-stage calculation method to convert 
the three-phase a, b, c measurements to the a and B two-phase model according to (2) and (3). And reverse this 
conversion, as shown in (2) and (3) [14]. 

Vb 


Cait 
lial = Filo = —=]| I: | 3) 


The real power P and the apparent power Q are calculated as shown in (4), (5) 


al=lve -vallis| @ 


AP=P+P (5) 


(2) 


where 


where AP is represents the average components of real power, P is AC component of P, P is DC component. 
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The output current in (7) represents the three-phase compensating current and it is an obligatory model of the 
three-phase PWM inverter, the above equation represents the application of the inverse Clark's transforms [15]. 


Figure 2. Block diagram of the instantaneous reactive power theory [5] 


3. HYSTERESIS BAND CURRENT CONTROLLER HBCC 

The hysteresis band current control technique is employed to control the compensated currents and 
select the switching signals for suitable VSI inverter gates. It is also because of the high accuracy, simple 
operation and independence of parameter variation [16]. The calculations of this model depend on the value of 
the error function in the current, which represents the difference between the source current and the 
compensated current injected into the network in (8) [17]. 


e(t) = Iref (t) — linj (t) (8) 


Figure 3 shows the working principle of this model on the comparison between the error function and 
the hysteresis band, as the difference between the reference current and the compensated current is more than 
the hysteresis band. The lower switching is turned on and the upper switching of the inverter is turned off, as 
a result, the current begins to fade away [18]. When the difference between the reference current and the 
compensated current is greater than the lower limit of the hysteresis band, the upper switching is turned on and 
the lower switching of the compensator is turned off, as a result, the current returns to the hysteresis band. The 
function (9) shows the conversion performance [19]. 


< f if Iinj > Iref + HB 


1 if Iinj < Iref — HB (9) 
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Figure 3. Hysteresis compartar 


4. DC LINK VOLTAGE CONTROL 

The P-type controller is used in the SAPF active power filter to regulate the DC-link voltage and the 
Lévy flight distribution algorithm is used to optimize the energy storage of the DC voltage in the capacitor 
based on the processing of the change in DC voltage E(t) to adjust the parameters of the controller Kp and Ki 
of in order control to progress the dynamic performance and obtain a stable state of the system, as shown in 
the Figure 4 [20]. 


C(s) = Kp += (10) 


where Kp represents the relative gain and provides the best dynamic performance, and Ki represents the 
integration gain that erases the error in the steady state. Use the traditional method the coefficients of the Kp 
and Ki controllers are calculated by comparing the system's total transfer function with the second-order 
generalized transfer [21]. 


Figure 4. Voltage control loop 


In the private accounts of the shunt system active power filter pi control is necessary to set the Piosses used in 
pq accounts as shown in (11) [22]. 


Plosses = KpAVdc + Ki f AVdc dt (11) 


where P loss represents the average loss that occurs inside the VSI transformer carried from the DC-bus link 
and is intended to give a good response and a suitable compensation where the real DC bus voltage (Vdc) is 
compared to the reference (Vdc*) and it is managed by the PI controller [23]. In this article, the target function 
is the error resulting from comparing the source Vdc and the voltage via DC-link. The control PI represents 
the flight position of the Lévy flight. The goal is to find the values of the constants Kp, Ki, and the error rate 
must be reduced [24]. 

In the mathematical model of the simulation environment [25], it is represented as a network 
environment of wireless sensors in its area. The algorithm's computational mechanisms start from calculating 
the Euclidean distance ED between each two adjacent sensor points and then determines if the sensor node is 
by the algorithm based on the calculated distance ED or moving it in another position. A new LFS model will 
be used elsewhere for a lower sensor node near the search spaces to reduce the occurrence of interference 
between the sensor nodes. Mathematically, the Lévy flight algorithm calculates the Euclidean distance ED 
between two adjacent sites, one being Xi and the other Xj using (12) 


ED(Xi; Xj) = y (xi — xj)? + Qi — yj)” (12) 
where (xi, yi) is the position coordinate of Xi and (xJ, yJ) is the XJ position coordinate. If the resultant distance 


is less than the threshold (this means there is a problem in deploying the agents in the search space), then the 
algorithm starts its mechanism by adjusting the positions of these agents using (13), 
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Xj(t + 1) = L'evy Flight (Xj (t); XleaderLB; UB) (13) 


The t is an index for the number of iterations. L evy_Flight is a function that performs the work of 
Lévy flights in terms of the step length and direction. Xj(t+1) illustrates the pseudo-code of this function. UB 
and LB are the highest and lowest values in the 2D dimensions of the search space. Respectively, the X Leader 
is the position of agent that has the lowest number of neighbors and will be used as the LF direction. Equation 
(14) moves the agent XJ towards the agent’s position that has the lowest number of neighbors. 


Xj(t + 1) = LB + (UB — LB)rand(n) (14) 


rand(n) is a function used to produce random numbersR in uniform distribution [0, 1]. In (15), offer 
more opportunities to discover solutions for unvisited sites in the search space and increases the exploration 
stage of the proposed algorithm. This equation updates the position of Xj to a new region in the search space 
where there are no other agents. 


R = rand(n).CSV = 0.5 (15) 


In each update for the Xj position, CSV is the comparison scalar value with R. The method examines 
the value of R in (15). At each iteration to update the sensor node XJ location (16). If the value is smaller than 
CSV in (14). If this is not the case, use (15) to give the algorithm additional chances to find the search space. 
The potential to explore and improve the algorithm's performance will be enhanced by varying the algorithm's 
solutions. The suggested approach uses to update Xi location. In (16) calculates Xi's new position, whereas 
(17) calculates Xi's ultimate position. 


Xi(t + 1) = TP + a1 x TFNeighbours + rand(n) x a2 x ((TP + a3XLeader)/2— Xi(t)) (16) 


XNewi (t + 1) = Levypigne(Xi(t + 1).TP.LB.UB) (17) 


TP is the position (solution) that obtained the greatest fitness value of the objective function, which is 
termed the target position, and Lévy flight is a function illustrated in (13). The random numbers 1, 2, and 3 are 
arranged in such a way that 0, 1, 2, 3 equals 10. TF neighbours surrounding Xi (t) have a total goal fitness of 
TF neighbours. 


TF NEIGHBOURS = ae 


(18) 
where Xk denotes Xi national's position, K denotes the neighbor's index, NN is the total number of Xi national's 
and D(K) denotes the fitness score for each vicinity. 


a1(V—Min(V)) 
Max(V)-—Min(V) 


D(k) = +02 (19) 


where, 


_ Fitness (Xj(t)) 
~~ Fitness (Xi(t))’ 


and 0 <@1, 02 < 1 (20) 


The suggested LFD method is depicted with the above equation. Each iteration of the LFD algorithm 
will repeat all of the previously stated equations. Assume that t represents the total number of iterations and n 
represents number of agents 


5. SIMULATION AND RESULT 

The simulation results for the designed model show an excellent compensation for the source 
harmonic current and an ideal compensation for SAPF. The model suggested had been implemented in a 
MATLAB simulation model, and a validation of the effectiveness of the designed PI controller in maintaining 
DC voltage across both ends of the capacitor with balanced load conditions. The specification parameters used 
in designing the model are listed in Table 1. The objective function is reduced by optimizing the PI parameters, 
the DC-link voltage control model is designed and the objective function is minimized by using optimization 
techniques, in this model, the Lévy flight determination algorithm was used, it is vital to set the search space 
to improve the parameters of the PI console, which are Kp, Ki, and get the best results. The search space can 
be specified for parameter 0.1< Kp<20, the search space for parameter Ki 10< Ki<200, the number of agents 
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10, and the number of iterations 10, we have run different simulations with the number of iterations to get an 
optimal value of the parameters of. Source stream analysis by fast Fourier series obtained from simulation in 
MATLAB. Simulation results for THD are included in Table 2 for the different cases before and after using 
the filter and when using the optimization method. The source current wave before filtering is shown in Figure 
5, and the load current before the filtering process in Figure 6 shows the extent of harmonic distortion. The 
THD of the source current before using the filter is 24.44% as shown in Figure 7 spectrum analysis of the 
source current. When using SAPF; While the compensator currents shown in Figure 8 that are pumped. The 
change in the wave form of the source current is obvious after filtering compared to the wave form of the load 
and filter current as shown in Figure 9. When using a conventional PI system, the THD decreases to 4.62% as 
shown in Figure 10. Also optimizing the operation of the PI control system, using a Lévy flight distribution 
algorithm graph, the THD of the source current is reduced to 1.56% as shown in Figure 11, and its waveform 
is as close as possible to the pure sine wave compared to its distortion. 


Table 1. Specifications of model design parameter 


Specifications Parameters Values 
Supply phase voltage (RMS) Vs 220 V, 50 Hz 
Line inductor Teac 3 mH 
Nonlinear load Diode bridge 63 A 
DC-link voltage Vac 620 voltage 
DC inductor Lee 1.6 mH 
DC capacitor Cac-equ 1.0 mF 
DC load: constant power resistive load Ry 40 Q 


Table 2. Simulation results of model 
Control technique Before using SAPF After using SAPF (Conventional PI) Using PI controller with (LFA) 


THD 24.44% 4.62% 1.56% 
Source Current 
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Figure 5. Source current without filter 
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Figure 6. Load current 


FFT analysis 
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Figure 7. Source current spectrum without filter 
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Figure 8. Filter current 
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Figure 9. Source current after using SAPF 


FFT analysis 
Fundamental (50Hz) = 9.431 , THD= 4.62% 
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Figure 10. Source current spectrum with filter (PI conventional) 


FFT analysis 
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Figure 11. Source current spectrum with filter (Lévy flights algorithm) 


6. CONCLUSION 

In this research, a group of techniques were used that showed their effects use and obtaining the best 
results in controlling the SAPF device. The pq theory was used in the current reference generation techniques 
for the speculative current calculations. The hysteresis band current controller HBCC theory was used to 
control the transformer device (VSI) within the PWM model. The Lévy flights distribution algorithm is 
proposed to improve the performance of the DC-link PI control, based on the techniques used in this model to 
determine the compensated reference current. The proposed model shows excellent performance in the SAPF 
device. These obtained results showed an increase in the power factor to 98.5 and that the voltage difference 
across the two ends of the DC. The capacitor is almost constant with a small ripple and that the quality of the 
source current is in an almost pure sinusoidal waveform after filtering, that the voltage waveform is in the same 
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phase, and that the voltage waveform is in the same phase, the THD was reduced in the current model from 
24.44% to 4.62% in the case of using the filter when improving the performance using the Lévy flight 
distribution algorithm, the THD was reduced to 1.56%, as the use of this method can be considered as one of 
the modern and important contributions that are used for the first time in this field and did not address it one 
of the researchers previously because the innovation of this algorithm is recent. 
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